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Abstract. The analysis of measurements of accelerated observers in Minkowski spacetime has led to the development 
of nonlocal special relativity theory. Inertia and gravitation are intimately connected in accordance with the principle of 
equivalence. We therefore seek a nonlocal generalization of the theory of gravitation such that in the new theory the field 
equations are integro-differential equations for the local gravitational field. We show that it is possible to develop a nonlocal 
generalization of Einstein's theory of gravitation via the introduction of a scalar "constitutive" kernel in the teleparallel 
equivalent of general relativity. The resulting nonlocal theory is essentially equivalent to Einstein's theory plus "dark matter". 
That is, nonlocality simulates dark matter by introducing a new source term into general relativity. In the linear approximation 
for the nonlocal modification of Newtonian gravity, we recover the theoretical basis for the phenomenological Tohline-Kuhn 
modified gravity approach to the explanation of the astrophysical evidence for dark matter. 
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INTRODUCTION 



A crucial observational aspect of gravitation is the universality of free fall within the classical domain of physics. 
In the context of Newtonian theory of gravitation, this experimental result implies the principle of equivalence of 
inertial and gravitational masses. The deep connection between inertia and gravity is the origin of Einstein's principle 
of equivalence, which establishes a definite local relationship between acceleration and gravitation |1]. Therefore, 
I following Einstein, to develop a theory of the gravitational field, our starting point must naturally be the physics of 

(~ — ■ accelerated observers in Minkowski spacetime. 

CO I Imagine a global background inertial reference system with Cartesian coordinates = (cf ,x). Henceforth we use 

units such that c = 1, unless specified otherwise. Consider an accelerated observer following a world line x'' (t) and an 
orthonormal tetrad frame A'^^aj caiTied by the observer along the path such that A'^(o) = dx^ /dz is the unit timelike 
vector tangent to the path as well as the four- velocity vector of the observer. Here T is the proper time along the world 
line. The motion of the local frame of the observer is given by 

■ where i^[a)(P) = ~'^{fi){a) is the antisymmetric acceleration tensor of the observer In analogy with electrodynamics, 

<p{a){p) has components (a, co), so that 0(a)(/3) can be expressed in terms of the translational acceleration of the observer 
(a) and the rate of rotation of the spatial frame with respect to a nonrotating (i.e., Fermi- Walker transported) frame 
(co). Out of these spacetime-invariant quantities, one can construct acceleration scales that characterize the variation 
of the state of the observer Once the kinematic aspects of the motion of an arbitrary observer in Minkowski spacetime 
are established, we turn to the problem of physics in accelerated systems, namely, what accelerated observers actually 
measure. 

In the standard theory of special relativity, Lorentz invariance is extended to accelerated observers in a pointwise 
manner via the hypothesis of locality [2-5]. That is, an accelerated observer is assumed to be pointwise inertial. 
This assumption originates from Newtonian mechanics, where the state of a point particle is characterized only by its 
position and velocity. For instance, the locality principle implies that the proper time t is in fact the time as determined 
by the accelerated observer ("clock hypothesis"). This approach is justified for phenomena involving classical point 
particles and rays of radiation, since in these regimes instantaneous measurements are in principle possible. It generally 
breaks down, however, for basic field measurements, as these cannot be performed instantaneously. Indeed, according 
to Bohr and Rosenfeld, the electric field E(f,x) and the magnetic field B(f,x) occur in Maxwell's equations as 
idealizations; only the spacetime averages of these fields have immediate physical significance l^lltl- While this 
observation appears innocuous for inertial observers, it acquires the status of a fundamental principle for accelerated 



observers due to the existence of invariant acceleration scales. The Bohr-Rosenfeld principle implies that in general 
field measurements cannot be instantaneous and hence must depend on the past world line of the observer fl- 

Let be a basic field in the global background inertial frame and let \jf = Ay/ he the field measured instantaneously 
by the infinite set of hypothetical momentarily comoving inertial observers along the world line of an accelerated 
observer Here A is a matrix representation of the Lorentz group. Suppose that 4'(t) is the field actually measured by 
the accelerated observer We are interested in the relationship between ^ and The most general relation between 
and 1// that is consistent with causality and linearity can be expressed as 

'I'(t) = v/(t) + m(t-to) f K{x,x')\if{T')dx' . (2) 

J To 

Here u{t) is the unit step function such that u{t) = for f < and u{t) = 1 for f > 0, Tq is the instant at which 
acceleration is turned on and the kernel K vanishes in the absence of acceleration. Equation (2) is a Volterra integral 
relation [9, -1 1,1 and represents the idea that what is measured at proper time T by the accelerated observer is the 
field V'^(t) measured by the instantaneously comoving inertial observer at T together with a certain average over the 
observer's past world line that constitutes the linear memory of past acceleration. 

For sufficiently low accelerations, the observer can be considered almost inertial and hence the nonlocal part of 
equation (2) constitutes only a small perturbation in this case. In general, if the intrinsic scale of the phenomenon 
under observation is very much smaller than the scale of variation of the state of the observer, then the deviation from 
locality is expected to be negligible. This is indeed the case for most Earth-based observations as c^/ja©! ~ 1 light 
year and c/jo)®! w 28 astronomical units. 

To determine the nonlocal contribution to equation (2), we need a physical principle that could form the basis for 
the calculation of the kernel. The general phenomenon of spin-rotation coupling furnishes the main idea: No observer 
can stay completely at rest with a basic radiation field. This is a generalization for all observers of a well-known result 
of Lorentz invariance for inertial observers. It follows from a detailed analysis that in general the kernel is directly 
proportional to acceleration and is given by lfT2l - [l6ll 

^(T,T')=fe(T') = -^A-'(T'). (3) 

For instance, this is the appropriate kernel for the Dirac field as well as for the electromagnetic gauge field {x). The 
special case of the electromagnetic field F^y{x) requires a separate treatment fl7*]. 

A nonlocal theory of accelerated observers, first proposed in 1993 [8], has been further developed in its essential 
aspects and the main outline of a complete nonlocal special relativity theory has been presented in 2008 [18]. 
Some of the observational consequences of the new theory have been worked out, especially regarding spin-rotation 
coupling lfT9ll . In this connection, Bohr's correspondence principle has been employed to show that the nonlocal theory 
is in better agreement with quantum mechanics than the standard local special relativity theory |15]. Moreover, the 
nonlocal theory forbids the existence of basic scalar or pseudoscalar fields. The nonlocal theory is consistent with all 
available observational data and its novel consequences should be subjected to direct experimental tests lfl9ll . 

It is important to emphasize that nonlocal special relativity actually involves local fields ^(x) in Minkowski 
spacetime that satisfy integro-differential field equations. It turns out that it is not possible to develop a consistent field 
theory of bilocal fields W{x,x'). Therefore, the fields are local, but the acceleration-induced nonlocality that originates 
from the memory of past observer acceleration appears in the circumstance that the local fields satisfy nonlocal field 
equations |20]. 

Once a proper nonlocal theory of accelerated systems is on hand, it is natural to seek a nonlocal theory of gravitation 
due to the basic connection between inertia and gravity. At first sight, it might seem that one could employ Einstein's 
principle of equivalence and arrive at nonlocal general relativity in a straightforward manner. It turns out, however, 
that Einstein's beautiful approach does not work in this case due to the inherently local nature of Einstein's principle 
of equivalence. For instance, suppose that we start with Einstein's general linear approximation to general relativity 
theory as a classical spin-2 field in Minkowski spacetime. Nonlocal theory of accelerated systems can then be applied 
in this case to generate a linear nonlocal spin-2 field theory as determined by accelerated observers in Minkowski 
spacetime ll2lll22ll . It is however not clear how to implement in this case Einstein's heuristic pointwise identification 
of an observer in a gravitational field with a certain accelerated observer in Minkowski spacetime. Thus a direct attempt 
at a nonlocal generalization of Einstein's theory of gravitation appears to be futile. 

It is possible to arrive at general relativity (GR) from the standpoint of the gauge theories of gravitation. Indeed, 
the gauge approach to gravity naturally leads to spacetime theories with curvature and torsion. There is a spectrum 



of such theories such that at one end of the spectrum, one has GR based on a Riemannian spacetime manifold with 
only curvature and no torsion, while at the other end of the spectrum are spacetime theories with torsion and no 
curvature Of the latter, there is a unique one that is essentially equivalent to Einstein's general relativity: 

This is the teleparallel equivalent to general relativity (GR||). Teleparallelism has a long history; its application to 
gravitational physics has been considered by many authors starting with Einstein in 1928 ll28ll . GR|| is a translational 
gauge theory of gravity and has been recently studied by a number of authors ll29l - l46ll . Spacetime translations form 
an Abelian group; this fact led Friedrich Hehl to suggest that one should attempt a nonlocal GRy in analogy with 
the nonlocal electrodynamics of accelerated observers |47] in order to arrive indirectly at a nonlocal generalization of 
general relativity. This fruitful suggestion has been developed in recent papers Ii,48,-,50il . The resulting nonlocal theory 
implies that gravity is nonlocal even in the Newtonian regime. The nonlocally modified Newtonian gravitation appears 
to provide a natural explanation for the dark matter problem; that is, nonlocality simulates dark matter. 

The plan of this paper is as follows. GRy is briefly described in the following section. Its nonlocal generahzation 
is then presented and the corresponding field equations are worked out, for the sake of simplicity, in the linear 
approximation. The nonlocal modification of Newtonian gravity is discussed and its connection with the dark matter 
problem is emphasized. We conclude with a brief discussion of these results. 



GR|| 

The arena for the teleparallel equivalent of general relativity is the Weitzenbock spacetime. This gravitational theory 
can be characterized as a tetrad theory. The gravitational potentials are given by the four coframe 1 -forms d-" = ei"dx\ 
where e," (jc) is the tetrad field. Here Latin indices range over 0,1,2,3 and refer to holonomic spacetime indices, while 
Greek indices range over 0, 1,2,3 and refer to anholonomic tetrad indices. The latter are raised and lowered via the 
Minkowski metric tensor 7]„|3 = diag(l, — 1, — 1, — 1). The frame vectors ea = adj are dual to the coframe 1-forms, 
so that 

e,V^ = 5^«, e,Va = 5/. (4) 
The gravitational field strength is given by C" :— di}", hence 

C" = ^Cij"dx' Adx\ (5) 

where Cij" can be expressed as 

C,/-2a[,e^.]«. (6) 
The spacetime interval is given by c/s^ = gijdx' (g) dx-', where 

8ij^e,"e/r]c,^. (7) 

We define a geodesic in this spacetime connecting two fixed events P' and P to be such that the corresponding interval 
is an extremum; that is, 5 fpi ds = 0. This implies that 

d^x' ^ f ' 1 '^^'^ _ Q (-g-j 
ds'^ \ik] ds ds ' 

where the Christoffel symbols are given by 

= ^g'' {gij,k+gik,j-gjKi) ■ (9) 

In the Weitzenbock geometry, the frame field is globally teleparallel. This means that the connection 1-form, 
r«/^ = -r'^" = Ti^'f^dx', is such that 

Ra^ dTa^-TjATy^ - ^-R,ja^ dx' Adx' =0. (10) 
Therefore, one can choose a suitable Cartesian tetrad frame such that 

r"^ = o. (11) 



The situation here is similar to that of Minkowski spacetime in which the Riemannian curvature vanishes; then, it is 
possible to choose inertial Cartesian coordinates such that the connection coefficients all vanish globally. Henceforth, 
a star over the equality sign indicates that the corresponding relationship is only valid in the particular gauge in which 
equation (11) holds. 

While the frame field is globally parallel and curvature-free, the spacetime torsion is nonzero. In fact, :— 
is simply the gravitational field strength C" in the special gauge (11). To express the gravitational field equations in 
terms of C,j", it is useful to define the modified gravitational field strength 

G:,/:=k,/-C"[,;]+2q,«C,]/. (12) 



Furthermore, we have 



The GRi I field equations are then given by 



r^ij ±y__Sfrtj (13) 



d[iCj,f=0, (14) 

dj^ja-^a'^^a', (15) 
where T^' = \f—gJ'a\ Ta is the energy-momentum tensor and £«' is defined by 

£«'■ := -^e'a(C/ +Ca,^J{% . (16) 

The gravitational field equations in this framework have an interesting analogy with electromagnetic field equations 
due to the fact that both are gauge theories of Abelian groups. If we first take £«' in equation (15) to the right-hand 
side and then ignore for the moment the index "a" in equations (14) and (15), we see that these equations are just 
like Maxwell's equations in a medium. That is, C,j is analogous to the electromagnetic field tensor, so that C,/ has 
components (E,B) and, similarly, has components (D,H), so that £' in effect has the interpretation of current 
due to the field itself, since the gravitational field equations are nonlinear Indeed, £«' is the energy-momentum tensor 
density of the gravitational field in GR||. In this interpretation, the tetrad field e,", when we ignore "a", is the analog 
of the electromagnetic vector potential (i.e., gauge field), except that there are in fact four such potentials here, each 
for a = 0, 1,2,3, as there are four spacetime translations. Let us recall here that GR|| is the gauge theory of the four- 
parameter group of translations in spacetime, whereas electromagnetism is the gauge theory of the one-parameter C/ ( 1 ) 
group. Just as in electrodynamics, the homogeneous field equation ( 14) is automatically satisfied as a consequence of 
equation (6), while equation (13) is the analog of the constitutive relation for the medium. 

It is important to note that the GRy field equations can also be derived from the action principle 

5 y (£g + X:m)^^'*x = 0, (17) 

where £g and are the gravitational and matter Lagrangian densities, respectively. Moreover, in analogy with 
electrodynamics, £g is quadratic in the gravitational field strength, so that 

'^g = ^<^'i"^^^ = -^^'^«<^'7" (18) 



and 5Z^/5e" = 7a . Further aspects of the GR|| field equations can be found, for instance, in ll28ll and ll49ll . 
The proof that in the classical domain this theory is indeed equivalent to GR is given, for instance, in ll28ll . 



NONLOCAL GR|| 

We recall that in the phenomenological electrodynamics of media, it is possible for the constitutive relation to be 
nonlocal; for instance, one can have 

D(0=E(f) + y' fit-t')E{t')dt' , (19) 



where / depends on the properties of the medium — see Eq. (58.3) of ifsill . Nonlocal electrodynamics of accelerated 
systems has a similar formulation ll47ll . This correspondence suggests a natural nonlocal extension of GR| | , namely, that 
we should maintain the field equations (14) and (15), but replace the local "constitutive" relation (13) with its nonlocal 
generalization in Weitzenbock spacetime. To this end, we must postulate the existence of a "constitutive" kernel that 
would be responsible for the nonlocal nature of the gravitational interaction. The simplest possibility would be a scalar 
kernel. Moreover, we must ensure that the nonlocal extension of equation (13), once it is fully expressed in terms of 
holonomic indices, would be consistent with general coordinate invariance. This can be done by means of the world 
function £2 and its covariant derivatives |52]. Let us assume that there exists a unique geodesic joining every pair of 
events x and x' in the spacetime region of interest; then, £2 is a smooth biscalar in (x,x') such that 2£2 is the square of 
the invariant geodesic interval between x and x' . It satisfies the basic partial differential equations 

2£2 = /''U,ii6=^'^n,n., (20) 

so that the vectors Q.a{x,x') dH/dx" and £2j(x,x') := dH/dx'' are tangents to the geodesic atx andx', respectively, 
and have the same length as the geodesic. Here indices a,b,c,... refer to x and i,j,k,... refer to x'. Differentiating 
equation (20), we find that Q.ai{x,x') = Q.ia{x,x') and these tend to —gai(x) in the limit where x and jc' coincide. Let us 
recall here that in Minkowski spacetime, 

aix,x') = ^ri,j{x"~x'){x'j~x') , (21) 

so that £2" = x" — x'", Q.' — x'' — x', £2„, — — rjai , ^ab — ^ab ™d £2,,- = r\ij . Further properties of the world function 
can be found in |52] and |49]. A possible nonlocal generalization of equation (13) is then 



^"'c{x) = ^^^^[e:"'^.(x)-|c/(x,x')£2«'£2^^£2,, 

xK{x,x')(Lij''{x')^J -g{x')d'^x'] . (22) 

Here the function U has been introduced to ensure that equation (22) is consistent with causality; that is, x' must be in 
the past of x. Thus t/ = if x' is in the future of x; otherwise, U ~\. 

It remains to discuss the nature of our nonlocal scalar kernel K{x,x'). It could depend upon £2 as well as spacetime 
scalars formed from its covariant derivatives such as, for instance, £2"ea"(x) and £2'e,"(x'). Moreover, K could depend 
on the invariants constructed from the gravitational field strength Cjjk- As a third rank tensor that is antisymmetric in 
two indices, Cijk has three irreducible pieces: a tensor, a vector and an axial vector. Out of these, one can build three 
hnearly independent quadratic invariants that turn out to be equivalent to the three Weitzenbock invariants 

CijkC'\ CkjiC'i', Q/C'k (23) 

at X and x'. Of course, one could also imagine various combinations of such possibilities. 

The substitution of equation (22) in equation (15) results in the main nonlocal gravitational field equation. The 
analysis of the physical consequences of the resulting nonlinear and nonlocal theory is a daunting task; therefore, we 
resort to the linear approximation in the rest of this paper. 



LINEAR APPROXIMATION 

The linear approximation in GR| | is defined by 

e'a = 5i,-w'a, ei"^5," + xif",, (24) 

where y/ja is assumed to be small and will be treated to first order, in which case the distinction between holonomic 
and anholonomic indices disappears. It follows that 

gij = riij + hij , hij = 2 Y^ij) , (25) 

so that the gravitational field strength is given by 

Q/ = 2w\j,], (26) 



(27) 



where x^r = 77,7 i^'A 

Let us introduce the trace-reversed potentials hn^ such that 

hik^hik-^riikh, h^riijh'' =2x1/. (28) 

Then, 

dj^j, = - ^nhk, + + - i iikih'' ji , (29) 

which is simply the symmetric Einstein tensor G,i in the linear approximation. We now need to substitute the linearized 
form of the nonlocal "constitutive" relation (22) in the linearized version of equation (15). Therefore, the nonlocal field 
equation can be written in this approximation as 

Gikix) + r^"'> I ^^^e,,,iy)d^y = 87tGUx) , (30) 

where K{x,x') := U{x,x')K{x^x') and diT^' = 0. Let us note that in the absence of the nonlocal kernel, we recover the 
linear approximation to general relativity. In general, Tij is not symmetric; therefore, in equation (30) we have 16 field 
equations for the 16 components of . 



NEWTONIAN LIMIT 

The Newtonian limit in the present framework corresponds to the assumption that 



1 24> 

V^OO = ^11=^22 = ^33 = — * , W= r ' f^^l) 



c- 



where <I> is the Newtonian gravitational potential and the other (off-diagonal) components of yz/y are neglected in this 
limit. This is related to the fact that in the Newtonian regime of GR, gij = Tj,y + /i,^ , where hij = 2c^^<P diag (1,1,1,1). 
We therefore have c^Coao ~ —^d^'P, where A, B, C, ... here indicate spatial indices that range over 1, 2 and 3. Moreover, 
c^Goo = 2V^<P as in GR and Too = pc^, so that equation (30) reduces to 

V^*W + E/^^<'^^-4.Gp(x), ,32, 

where 

K{x,y) = 5ix°-y'')k{x,y), (33) 

since all retardation effects are neglected in the Newtonian limit (c °°). 

The kernel A:(x,y) in this nonlocal modification of Poisson's equation of Newtonian gravity can depend on x — y 
as well as on the Weitzenbock invariants at x and y. The latter have been calculated in this case |5Qil and the result is 
essentially (VO)^. Thus the kernel can be written as 

Mx,y)=fc'(x-y)+fc"(^x-y;^) . (34) 

The second term in equation (34) involves a nonlinear modification of Newtonian gravity. To simplify matters, we 
concentrate on the first term in what follows. The resulting linear nonlocal modification of Poisson's equation is then 



V24>(x) + k'(x- y)Vl^iy)d^y = 47rGp (x) . 



(35) 



Using the standard Liouville-Neumann method of successive substitutions, equation (35) can be expressed as 



V^<I> = 4;rG(p + pD), (36) 

where 

Pd{x)^ [q{x-y)p{y)d'y. (37) 



Here q' is a universal function that is independent of the nature of the source and is the reciprocal kernel to k' 
Thus nonlocality can appear as an extra matter source ("dark matter") in the standard local theory. 

It is interestin g to see if the nonlocal theory can actually provide a natural solution to the dark matter problem in 
astrophysics ||53| - d5I1 . To this end, we consider a spiral galaxy of mass M and, for simplicity, let p(x) =Md{x) and 

be the density of dark matter, with corresponding asymptotic speed vq, that can account for the flat rotation curves of 
spiral galaxies. Then it follows from equation (37) that 



where X ~ GM/vq is a constant of order 1 kpc for spiral galaxies. The solution of the modified Poisson equation in 
this case is 

GM GM^ /|x|\ 
4>(x) = -^ + — in U . (40) 



It turns out that such a modification of Newtonian gravity was first phenomenologically proposed by Tohline 156|] and 
later developed by Kuhn and his coworkers IstI 15811 in order to account for the dark matter problem in spiral galaxies 
and galaxy clusters. In the Tohline-Kuhn scheme, A is a universal constant and hence M °^ v^; however, the empirical 
Tully-Fisher relation lf59l l60l] implies that M o= Vq is observationally preferred. The inclusion of nonlinearity in the 
modified Poisson equation via an appropriate ^" in equation (34) may solve this problem. 

The implications of modified Newtonian potential (40) for gravitational physics in the solar system are briefly 
studied in |50]. It is important to look for any observational evidence of nonlocal gravity in the solar system; however, 
at present such effects appear to be too small to be detectable ||50|] . 



DISCUSSION 

Starting from first principles, we have presented a framework for nonlocal gravity that has the potential to provide a 
natural theoretical explanation for the observational evidence regarding dark matter. Indeed, nonlocality appears to 
simulate dark matter by introducing in effect a new source term ("dark matter") in Einstein's theory of gravitation. 
In the Newtonian regime, the density of dark matter is essentially the convolution of matter density with a kernel 
that represents the novel nonlocal aspect of the gravitational interaction and must ultimately be determined from 
observation. A more detailed treatment of nonlocal gravity is contained in l,48. - .50.] . 
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